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(57) Abstract 

A method of applying scatter 
and attenuation correction to emis- 
sion tomography images of a region 
of interest of a patient under obser- 
vation comprises the steps of aligning 
a three-dimensional computer model 
representing the density distribution 
within the region of interest with the 
emission tomography images (56, 58) 
and applying scatter and attenuation 
collection to the emission tomogra- 
phy images using the aligned computer 
model as a guide (60). 
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APPLICATION OF SCATTER AND ATTENUATION 
CORRECTION TO EMISSION TOMOGRAPHY 
IMAGES USING INFERRED ANATOMY FROM ATLAS 

Technical Field 

The present invention relates to emission tomography and in particular to a 
method and apparatus for applying scatter and attenuation correction to emission tomography 
images using anatomy inferred from an atlas. 

Background Art 

Single Photon Emission Computed Tomography (SPECT) and Positron 
Emission Tomography (PET) are nuclear medicine diagnostic imaging techniques used to 
measure the three-dimensional distribution of a radiopharmaceutical within the body. Brain 
SPECT and PET imaging techniques are primarily used to measure regional cerebral blood 
flow in a patient injected with a radiopharmaceutical to assist in the evaluation of stroke and 
the diagnosis of dementias such as Alzheimer's disease. 

Although SPECT and PET are useful imaging techniques, their poor 
quantitative accuracy has been an obstacle in the abiUty to achieve increased diagnostic 
reliability. Quantitative accuracy of brain SPECT and PET imaging has however been 
improved significantly through the application of scatter and attenuation correction to SPECT 
and PET brain images. To be sufficiently accurate, the application of scatter and attenuation 
correction to SPECT and PET brain images must be guided by the distribution of density 
within the head. Unfortunately, the distribution of density within the head cannot be obtained 
from SPECT and PET brain scans and therefore, separate measurements are required. 

Transmission imaging has been used to measure the distribution of density 
within the head to allow scatter and attenuation correction to be applied to SPECT and PET 
brain images. Unfortunately, the hardware necessary for making transmission measurements 
is complex, unreliable and requires extensive maintenance. Also, the need to make 
transmission imaging measurements in addition to the SPECT or PET brain images, increases 
the time required to complete the overall imaging procedure. SPECT and PET imaging 
procedures are themselves lengthy and require a patient to remain motionless to ensure 
accurate brain images. For sick and elderly patients, this is a difficult task. Adding to the 
length of the imaging procedure increases the likelihood that patients will not remain 
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motionless. Movement of a patient during the transmission imaging procedure results in 
inaccurate measurements of the distribution of density within the head. This of course 
provides an inaccurate guide for the application of scatter and attenuation correction to 
SPECT and PET brain images. Accordingly, improved methods to increase the diagnostic 
reliability of emission tomography images are desired. 

It is therefore an object of the present invention to provide a novel method and 
apparatus for applying scatter and attenuation correction to emission tomography images. 

Disclosure Of The Invention 

Broadly stated, the present invention provides a method and apparatus for 
applying scatter and attenuation correction to emission tomography images which estimates 
or "infers" the distribution of density within a region of interest of a patient under observation 
using a three-dimensional computer model of the region of interest. It has been found that 
scatter and attenuation correction guided by a computer model of the region of interest xmder 
observation produces results similar to those when using transmission images of the actual 
region of interest as the guide to the application of scatter and attenuation correction. 

According to one aspect of the present invention there is provided a method of 
applying scatter and attenuation correction to emission tomography images of a region of 
interest of a subject under observation comprising the steps of: 

aligning a three-dimensional computer model representing the density 
distribution within said region of interest with said emission tomography images; and 

applying scatter and attenuation correction to said emission tomography 
images using said aligned computer model as a guide. 

In a preferred embodiment, the computer model is in the form of a two- 
component atlas. During the aligning step, a functional component of the atlas is firstly 
aligned with the emission tomography images to generate a set of spatial transformation 
parameters. Following this, an anatomical component of the atlas is ahgned with the 
emission tomography images using the set of spatial transformation parameters. 

The atlas may be selected from a database of atlases based on degree of 
registration with the emission tomography images. Alternatively, multiple atlases maybe 
combined to yield a resultant atlas which better registers with the emission tomography 
images. 
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According to another aspect of the present invention there is provided in an 
emission tomography imaging system where emission tomography images of a region of 
interest of a subject are taken for analysis and are corrected for scatter and attenuation, the 
improvement comprising: 

using a three-dimensional computer model approximating the density 
distribution within the region of interest as a guide to the application of scatter and 
attenuation correction. 

According to still yet another aspect of the present invention there is provided 
an emission tomography image processing system comprising: 

memory storing emission tomography images of a region of interest of a 

subject; 

said memory also storing at least one three-dimensional computer model of 
said region of interest, said computer model representing the density distribution within said 
region of interest; and 

a processor for registering said computer model with said emission 
tomography images and for applying scatter and attenuation correction to said emission 
tomography images using said registered computer model as a guide. 

According to still yet another aspect of the present invention there is provided 
an emission tomography imaging system comprising: 

means for taking emission tomography images of a region of interest of a 
subject to form a three-dimensional image of said region of interest; 

memory to store said emission tomography images, said memory also storing 
at least one three-dimensional computer model of said region of interest, said computer model 
representing the density distribution within said region of interest; and 

a processor for aligning said computer model with said emission tomography 
images and for applying scatter and attenuation correction to said emission tomography 
images using said aligned computer model as a guide. 

According to still yet another aspect of the present invention there is provided 
a computer readable medium including computer program code for applying scatter and 
attenuation correction to emission tomography images of a region of interest of a subject, said 
computer readable medium including: 
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computer program code for aligning a three-dimensional computer model 
representing the density distribution within said region of interest with said emission 
tomography images; and 

computer program code for applying scatter and attenuation correction to said 
emission tomography images using said aligned computer model as a guide. 

The present invention provides advantages in that by using a three- 
dimensional computer model of the region of interest of a subject under observation that 
approximates its density as a guide for the application of scatter and attenuation correction to 
emission tomography images, the need for transmission imaging is obviated. Therefore, in 
the case of SPECT and PET imaging, the imaging procedures do not need to be lengthened. 
Also, since the distribution of density in the region of interest under observation is 
approximated by a three-dimensional computer model, additional hardware is not required to 
create the guide for the application of scatter and attenuation correction. This makes the 
present invention significantly less expensive and more flexible than transmission imaging 
systems. In addition, since a three-dimensional computer model of the region of interest 
under observation is used as the guide for the application of scatter and attenuation correction, 
scatter and attenuation correction can be applied retrospectively to existing databases which 
include significant numbers of emission tomography images for which no transmission 
imaging measurements were acquired. 

Brief Description Of The Drawings 

An embodiment of the present invention will now be described more fully 
with reference to the accompanying drawings in which: 

Figure la shows a two-dimensional emission tomography brain image without 
the application of scatter and attenuation correction; 

Figure lb shows the two-dimensional image of Figure la with the application 
of scatter and attenuation correction; 

Figure Ic is a two-dimensional transmission brain image; 

Figures Id and le are two-dimensional emission tomography brain images 
with appurtenant anatomy derived from transmission images applied thereto; 
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Figure 2 is a block diagram showing a method for applying scatter and 
attenuation correction to emission tomography images in accordance with the present 
invention; 

Figure 3a shows a two-dimensional emission tomography brain image of a 
head phantom with scatter and attenuation correction; 

Figure 3b is a two-dimensional transmission brain image; 

Figures 3c and 3d are two-dimensional emission tomography brain images 
with appurtenant anatomy derived from transmission images applied thereto; 

Figures 3e and 3f are two-dimensional emission tomography brain images 
with inferred anatomy derived from an atlas applied thereto; 

Figures 4 and 5 show two-dimensional brain images and a graph comparing 
uniform appurtenant anatomy, uniform inferred anatomy and brain contour; 

Figure 6 shows two-dimensional brain images with non-uniform appurtenant 
and non-uniform inferred anatomy applied thereto; 

Figures 7 and 8 show quantitative evaluation of head phantom SPECT 
reconstructions; 

Figure 9 shows two-dimensionl SPECT brain images; 

Figure 10 shows SPECT brain images superimposed onto transmission 
reconstruction and inferred anatomy; 

Figure 11 shows profiles taken through reconstructions of a SPECT brain 

image; 

Figure 12 is a graph showing a comparison of regional cerebral blood flow 
from reconstructed SPECT images guided by inferred anatomy and transmission scans; and 

Figure 1 3 is a graph showing the correlation between the reconstructed SPECT 
images of Figure 12. 

Best Mode for Carrying Out The Invention 

During emission tomography imaging such as for example SPECT and PET, a 
patient is injected with a radiopharmaceutical. Two-dimensional images or projections of a 
region of interest (ROI) of a patient are taken along an arbitrary axis. The projections are 
digitized and conveyed to a computer for storage in a three-dimensional array in computer 
memory to reconstruct a three-dimensional image of the region of interest. Once the three- 
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dimensional image is reconstructed, two-dimensional slices of the region of interest can be 
viewed along virtually any arbitrary axis using conventional software. Figure la shows a 
two-dimensional emission tomography brain image. 

During the imaging procedure, the emissions from the radiopharmaceutical are 
scattered and/or attenuated by different density tissue, air cavities and/or bones in the region 
of interest under observation. As can be seen in Figure la, scattering and attenuation of 
radiopharmaceuticals affects quantitative image quality. As a result, acquired emission 
tomography images are often unreliable. Applying scatter and attenuation correction to 
emission tomography images, using transmission images of the same region of interest taken 
during the same imaging procedure, is known but suffers from the disadvantages discussed 
previously. Figure lb shows the two-dimensional image of Figure la with scatter and 
attenuation correction applied using transmission images taken of the same anatomy. Figure 
Ic shows an example of a transmission image and Figures Id and le show appurtenant 
anatomy (sagittal and transverse respectively) derived from the transmission images apphed 
to emission tomography brain images. 

In the present invention, a three-dimensional computer model or "atlas" of the 
region of interest, which provides accurate density distribution of the region of interest, is 
stored in a database in computer memory and is used as a guide for the application of scatter 
and attenuation correction to emission tomography images. In the present embodiment, the 
atlas includes two components, namely a functional component simulating a SPECT OR PET 
scan of the region of interest and an anatomical component simulating a transmission scan of 
the region of interest. The atlas can be created from existing transmission images or x-ray CT 
scans of similar regions of interest from other patients and averaged to form a suitable 
computer model or atlas. Multiple models for each region of interest and models for different 
regions of interest can be stored in the database and accessed individually during scatter and 
attenuation correction of emission tomography images. Since a computer model of the region 
of interest is used, additional hardware and procedure time is not required to apply scatter and 
attenuation correction to emission tomography images. An example of the application of 
scatter and attenuation correction to emission tomography brain images in accordance with 
the present invention will now be described. 

Referring now to Figure 2, a block diagram illustrating the present method for 
applying scatter and attenuation correction to emission tomography brain images is shown. 
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Initially, emission tomography brain images of a patient (block 50) are acquired in a known 
manner. The acquired brain images are in the fomi of two-dimensional projections of the 
radiopharmaceutical distribution in the brain. 

Initially, a preliminary reconstruction of the acquired brain images is 
performed to digitize and assemble the two-dimensional projections into a three-dimensional 
array in computer memory to create a three-dimensional image of the brain (block 52). An 
atlas of the head (the region of interest in this case) is then downloaded from the database into 
a datafile (block 54) and the functional component of the head atlas is identified. Once 
identified, the functional component of the head atlas, in this embodiment the brain 
component, is copied to a three-dimensional array. Following this, an alignment procedure to 
register spatially the brain component of the head atlas with the three-dimensional brain 
image is performed (block 56). In the present embodiment, a simplex algorithm is used to 
register spatially, the brain component of the head atlas with the brain image as described in 
"Numerical Recipes In C" by Press et al, 2""* edition, New York, New York, Cambridge 
University Press, 1992. Those of skill in the art will however appreciate that other 
registration algoritlims can be used to align the brain component of the head atlas to the three- 
dimensional brain image. During this ahgnment procedure, a set of 3D spatial transformation 
parameters representing the three-dimensional transformations, including but not limited to 
rotation, shifting and scaling, that are necessary to register the brain component of the head 
atlas with the three-dimensional brain image, is calculated and is stored in the computer 
memory as a matrix. 

Once the set of 3D transformation parameters is calculated and stored, the 3D 
spatial transformation parameters in the set are applied to the anatomical component of the 
head atlas to register it with the three-dimensional brain image (block 58). With the 
anatomical component of the head atlas aligned with the three-dimensional brain image, the 
atlas can be used as a density distribution guide to the application of scatter and attenuation 
correction to the three-dimensional brain image. 

With an accurate density distribution guide established, scatter correction is 
applied to the three-dimensional brain image followed by attenuation correction in a known 
manner (block 60). Once scatter and attenuation correction has been applied to the three- 
dimensional brain image, the brain image is reconstructed into a three-dimensional array in 
computer memory to complete the image correction process (block 62). The corrected three- 
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dimensional brain image can then be analyzed (block 64). Appendix A includes pseudocode 
representing the above-identified process. 

Figures 3a to 3f show a comparison of two-dimensional brain images of a head 
phantom with anatomy derived from transmission images (appurtenant anatomy) and 
anatomy derived from a head atlas (inferred anatomy) applied. 

The present method and apparatus was tested using an anthropomorphic head 
phantom modeling soft tissue, hard tissue and air cavities within a skull and including a two- 
compartment brain reservoir. The two compartments of the reservoir were separately filled 
with two water solutions of Tc-99m, having a specific activity ratio of 4:L Fan-beam SPECT 
was acquired followed by a Tc-99m transmission scan sixty hours later. The reconstructed 
transmission image is referred to as appurtenant anatomy. 

Five scatter and attenuation correction schemes were evaluated based on non- 
uniform appurtenant anatomy, uniform appurtenant anatomy, non-uniform inferred anatomy, 
uniform inferred anatomy, and uniform brain contour. For uniform scatter and attenuation 
correction, the inferred and appurtenant anatomies were segmented and assigned uniform 
attenuation coefficients of soft tissue (0.1 5cm"^ for Tc-99m). A SPECT reconstruction of the 
head phantom was similarly processed to facilitate scatter and attenuation correction guided 
by brain contour. Scatter correction was based on a non-stationary deconvolution scatter 
subtraction as described in the paper authored by Stodilka et al entitled "The Relative 
Contributions of Scatter and Attenuation Correction Toward Brain SPECT Quantification", 
Phys Med Biol, 1998. Attenuation correction/reconstruction was subsequently performed by 
ordered-subsets expectation-maximization as set out in the paper authored by Hudson HM et 
al entitled "Accelerated Image Reconstruction using Ordered Subsets of Projection Data", 
IEEE Trans Med Imaging 13 601-609, 1994. Although particular examples of scatter 
correction, attenuation correction and reconstruction algorithms are described, those of skill 
in the art will appreciate that other scatter correction, attenuation correction and 
reconstruction methods can be used. 

Uniform scatter and attenuation correction should be guided by the contour of 
the attenuating medium if an acceptable level of accuracy for objective diagnosis is to be 
expected. Figures 4 and 5 show that uniform appurtenant anatomy is better approximated by 
uniform inferred anatomy than by brain contour. As can be seen in Figure 6, non-uniform 
appurtenant anatomy and non-uniform inferred anatomy are similar. 
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Head phantom SPECT reconstructions were quantitatively evaluated by four 
metrics, namely bias, uniformity, contrast-recovery, and relative quantification (see Figures 7 
and 8). As will be appreciated, scatter and attenuation correction guided by inferred anatomy 
provides quantitative accuracy that is distinctively superior to scatter and attenuation 
correction guided by brain contour. 

Application of scatter and attenuation correction to emission tomography 
images of anatomy other than the head can also be performed. For example, inferred 
anatomy can also be used to apply scatter and attenuation correction to cardiac images. 
During construction of the cardiac atlas, an anatomical component of the cardiac atlas 
representing the anatomical features of the thorax is created and includes: 

soft-tissues such as the heart, liver, muscle, and fat; 

very low-density soft-tissues such as the lungs; and 

high-density tissues such as bone and cartilage in the ribs and spine. 

A functional component of the cardiac atlas is also created and simulates a 
SPECT or PET cardiac image. Appropriate data to construct the cardiac atlas can be obtained 
in a variety of ways including, but not limited to, imaging a phantom or human subject by X- 
ray CT, MRl, or gamma-camera transmission computed tomography; or computer simulation. 
Also, the cardiac atlas can be constructed by amalgamating a plurality of patient scans. 

The procedure of using inferred cardiac anatomy to apply scatter and 
attenuation correction to cardiac images remains the same as with brain imaging described 
above. First, a preliminary reconstruction of the patient's cardiac SPECT data is performed. 
The functional component of the cardiac atlas is then registered to the preliminary 
reconstruction. The registration includes a spatial transformation that may include shifting, 
rotation, scaling, and/or non-linear operations such as warping. The registration procedure 
calculates a matrix representing the spatial transformation that maps atlas space into patient- 
specific space. Once the matrix has been calculated, the matrix is applied to the anatomical 
component of the cardiac atlas, thus inferring anatomy in the chest. 

Although the above methods describe the use of a single generic atlas of the 
anatomy under observation, those of skill in the art will appreciate that custom atlases can be 
developed and stored in the database. For example, disease specific atlases such as an 
Alzheimer's disease atlas or a stroke atlas can be developed and used when correcting 
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emission tomography images of patients suffering from those diseases. The disease specific 
atlases may be tracer or lesion specific to allow for local concavities in radiopharmaceutical 
uptake. Patients with severe Alzheimer's or Pick's disease do not have normal cerebral blood 
flow. Areas with flow deficits can limit the accuracy of the registration of the emission 
tomography images with the functional component of the head atlas if the head atlas assumes 
normal blood flow and hence uniform radiopharmaceutcial uptake. Atlases can also be 
developed to take into account physical traits such as for example, an exceptionally large 
nasal sinus. During scatter and attenuation correction, the operator can select the appropriate 
atlas to use. Alternatively, the atlas can be selected automatically by computer software. In 
this case, the software, performs a preliminary reconstruction using each custom atlas and 
registers the atlas to the preliminary reconstruction to measure the accuracy of the 
registration. The atlas with the highest registration accuracy is then selected. If desired, the 
software can use fuzzy logic, theoretical calculations or other criteria to combine two or more 
atlases to yield a single resultant atlas, which provides a better degree of registration. 

EXAMPLE 

Example 1: Inferred Anatomy and Brain Imaging 

The following example is described for the purposes of illustration and is not 
intended to limit the scope of the present invention. 

Inferring Anatomy From A Head Atlas 

A head atlas was prepared as follows. A Zubal three-dimensional digitized 
MRI head phantom [Zubal et al 1994] was segmented to produce two data sets, namely a 
SPECT atlas simulating a SPECT scan of the phantom, and an anatomical atlas simulating a 
transmission scan of the phantom. The SPECT atlas consisted of voxels containing gray- 
matter and white-matter, to which ^^"Tc-HMPAO relative uptakes of 4 and 1 had been 
assigned respectively. The anatomical atlas consisted of voxels containing hard-tissue, soft- 
tissue and nasal sinus, to which the correspondent 140keV narrow-beam attenuation 
coefficients of 0.25, 0.15, and 0.075 cm"^ were assigned. 
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A patient SPECT scan was then reconstructed without scatter and attenuation 
correction. Facial activity was removed to yield a data set referred to as the preliminary 
patient reconstruction. The SPECT atlas was then registered to the preliminary patient 
reconstruction and the spatial transformation was recorded. This transformation was then 
applied to the anatomical atlas, thus inferring the location of the patient's soft-tissue, hard- 
tissue and air cavities (see Figure 9). 

A general-purpose radiological analysis program (Hermes, Nuclear 
Diagnostics, Stockholm, Sweden) was used to perform the unimodality registration. The 
large variation in head orientation necessitated that a manual registration be first performed. 
This was followed by an automated refinement. The cost function of the automated 
registration was defined as the sum of absolute count differences [Hoh et al 1993]. A global 
minimum was sought by a simplex search [Nelder and Mead 1965] within a parameter space 
consisting of rotating, shifting, and linear scaHng in x, y, and z directions [Holman etal 1991, 
Slomka et al 1995]. 

Sequential Transmission And Emission Imaging 

Ten dementia patients (5 females and 5 males, with a mean age of 64.3 years) 
were analyzed. For each patient, a transmission scan was first acquired. Patients became 
very relaxed during the quiet transmission scan, and were then injected with 740 MBq of 
^^"Tc-HMPAO. The SPECT procedure was started approximately 5 minutes post-injection. 
The SPECT system, which has transmission capabihties, is described in detail [Kemp et al 
1995]. It consists of a General Electric 400AC/T gamma camera (General Electric, 
Milwaukee, WI) with a 409.6 mm diameter circular field-of-view. Projections were acquired 
through a fan-beam collimator (Nuclear Fields, Des Plaines, IL) having a 600 mm focal 
length and 1.5 mm flat-to-flat hexagonal hole width. The transmission component includes a 
fi-ame mounted onto the camera's collimator that holds a tantalum-coUimated ^^""Tc line 
source along the focal Une of the fan-beam collimator. Collimation of both the line source 
and camera minimizes scatter. As a result, the transmission system effectively measures 
narrow-beam attenuation coefficients [Tsui et al 1989, Kemp et al 1995]. SPECT and 
transmission scans were acquired with a 20% energy window, centered on the ^^""Tc 
photopeak of 140 keV. The scans consisted of 128 projections, equally spaced over 360"*. 
Each circular projection was acquired into a 128 x 128 pixel square matrix (1 pixel = 3.2 
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mm). Both transmission and SPECT scans were 10 seconds per projection and count rates 
were approximately 70 and 1 .5 kcounts per second, respectively. All scans were corrected for 
uniformity using 100 million count flood images, and transmission scans were normalized to 
50 million count blank images. Radii of rotation varied among the patients; the smallest 
being 170 mm and the largest being 205 mm. Prior to reconstruction, all scans were rebinned 
to object-plane parallel-hole geometry via two-dimensional cubic interpolation. 

Scatter And Attenuation Correction And Reconstruction 

The SPECT data were reconstructed using a maximum-likelihood estimator 
with an unregularized 32-level ordered subset [Hudson and Larkin 1994] implementation of 
the expectation maximization algorithm [Shepp and Vardi 1982, Lange and Carson 1984] 
(OSEM). The four projections that were used per sub-iteration were equally spaced about 
360*". Attenuation was modelled in the matched projector/backprojector pair, and a scatter 
estimate [Stodilka et al 1998b] was added as an priori background following forward 
projection [Lange and Carson 1984, Bowsher et al 1996, Kadrmas et al 1998]. Both scatter 
and attenuation modelling incorporated the narrow-beam attenuation coefficients from 
transmission imaging or inferred anatomy. Detector response was not included. Four 
iterations of OSEM were used, following initialization with a uniform [Nunez and Llacer 
1990] support prior derived from transmission reconstruction or inferred anatomy. 
Reconstructions were then post-filtered [Nunez and Llacer 1990] using a three-dimensional 
Butterworth filter with an order of 8 and cutoff at 0.42 cm *. The transmission data were also 
reconstructed by the emission OSEM algorithm, following blank scan normalization and log- 
transformation. 

Line source coUimation, coupled with limitations in detector count rate 
capability and patient compliance resulted in less than ideal transmission statistics. To reduce 
the effects of transmission imaging noise propagation into the SPECT reconstruction [Xu et 
al 1991, Tung and Gullberg 1994], the transmission reconstmctions were segmented as 
follows. Soft-tissue in the reconstructed transmission volumes was forced to have uniform 
density. First, a large region of interest (ROI) was drawn around soft-tissue regions, and 
mean and variance estimates were calculated. Then, all voxels having count densities within 
±2 standard deviations of this mean were assigned to 0.15cm '. Thus, the transmission 
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reconstructions were characterized by noiseless soft-tissue, yet featured hard-tissue and air 
cavities. 

Template-Based Quantification 

Previous work has identified that a major confound to reproducible 
quantification, originates from manual and threshold-dependent placement of anatomical 
ROIs onto SPECT scans [Msaki et al 1998]. To reduce this subjective source of error, a 
normal template [Msaki et al 1998] onto which twelve bilateral volumetric ROIs are 
demarcated [Karbe et al 1994] was used. The ability to store the template and its ROIs 
ensures reproducibility of the analysis. This quantification procedure also introduces 
standardization to the analysis, which facilitates the exchange of data among different 
institutions [Evans et al 1988]. All reconstructed scans were registered to the normal 
template, herein referred to as "spatial normahzation". Prior to spatial normalization, voxels 
previously identified as facial activity were set to zero. Following superposition of the 
template ROIs onto each scan, the cortical rCBF for each ROI was normaUzed to cerebellar 
rCBF [Karbe et al 1994] and corrected for blood flow-dependent tracer reflux [Lassen et al 
1988]. Analysis of the absolute concentration of radiopharmaceutical was not performed 
since currently, absolute rCBF quantification is seldom used in SPECT [Bakker and Pauwels 
1997]. 

Quantitative Error Analysis 

A sample of inferred anatomy and transmission reconstruction is illustrated in 
Figure 10. In Figure 11, profiles through the SPECT reconstructions guided by transmission 
scans and inferred anatomy are presented. The profiles have not been normalized to 
cerebellar count density. 

The means and standard errors of regional cerebral blood flow from SPECT 
scans guided by inferred anatomy and transmission scans were compared ROI-by-ROI (see 
Figure 12). Statistical analysis was also performed ROI-by-ROI using repeated analysis of 
variance to determine where there were significant differences between the inferred anatomy 
versus transmission-guided reconstruction methods. Pooled-ROI and ROI-dependent 
correlation coefficients were calculated between inferred anatomy reconstruction and 
transmission-guided SPECT reconstructions. Figure 13 shows the correlation between the 
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two reconstruction and quantification methods by pooling all ROIs and patients together. 
The p value for significance was set to 0.05 for all tests. 

Inferred Anatomy Error Propagation Analysis 

Four sources of error were identified that contribute to discrepancies in ROI 

quantification: 

(1) inferred anatomy-guided scatter correction; 

(2) inferred anatomy-guided attenuation correction; 

(3) inferred anatomy-guided spatial normalization, which is equivalent to ROI 
misplacement; and 

(4) patient motion between transmission and SPECT scans. 

The first two are inherent limitations to the principle of inferred anatomy, 
whereas the last two represent artifactual exaggerations of errors in the context of this 
example. The first three sources of error, namely scatter, attenuation, and ROI misplacement 
were measured. 

The propagation of ROI quantification to differences between inferred 
anatomy scatter correction and transmission-guided scatter correction were analyzed by 
performing: 

(1) inferred anatomy-derived scatter correction; 

(2) transmission-derived attenuation correction; and 

(3) applying the spatial normalization calculated to be optimal for registering 
the transmission-derived SPECT reconstruction to the quantification template. 

A similar analysis for evaluating the effects of inferred anatomy-derived 
attenuation correction was carried out by performing: 

(1) transmission -derived scatter correction; 

(2) inferred-anatomy derived attenuation correction; and 

(3) applying the spatial normalization calculated to be optimal for registering 
the transmission-derived SPECT reconstruction to the quantification template. 
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The effects of ROI misplacement were quantified by perfomiing transmission- 
derived scatter and attenuation correction, and then applying the spatial transformation 
calculated to be optimal for registering the inferred anatomy-derived SPECT reconstruction 
to the quantification template. Thus, the full propagation analysis resulted in three 
reconstructions, each of v^hich was quantitatively compared (via the above-described 
template-based quantification procedure) with the "gold-standard" transmission-derived 
SPECT reconstruction and spatial normalization. This procedure was performed for the ten 
patients, and the results averaged. However, once the errors due to scatter, attenuation, and 
ROI misplacement were separately quantified, their totals were found to be less than the 
errors caused by the full inferred anatomy protocol. This additional source of error is termed 
"unaccountable" in Table 2, and is believed to be caused by patient motion. 

Qualitative Analysis 

A sample comparison of inferred anatomy and transmission reconstruction is 
shown in Figure 10 for mid-sagittal and cortical-axial slices. Good reproduction of soft-tissue 
and hard-tissue at the cortical level is noted. Some discrepancy is seen near the vertex; 
however, this region is seldom included in quantitative analysis. Discrepancies near the nasal 
sinus are most marked. Fortunately, these structures mostly involve low-density areas such 
as air, and to a much lesser degree, soft-tissues and cartilage, which do not scatter and 
attenuate photons as much as higher density structures. 

Mid-brain profiles, shown in Figure 11, compare a SPECT reconstruction 
guided by transmission imaging with the same SPECT reconstruction guided by inferred 
anatomy. The profile was taken along the longest axis of the head, which is most sensitive to 
mis-registered transmission maps following scatter and attenuation correction [Huang et al 
1979]. 

Quantitative Error Analysis 

Table 1 below shows the results of the repeated analysis of variance and 
correlation analysis comparing transmission reconstruction and quantitative and inferred 
anatomy-guided reconstruction and quantification. Left frontal and central sulcus ROIs 
showed the highest probability of a true difference, reaching statistical significance with 
p=0.001 and 0.002, respectively. These ROIs also had the highest correlation coefficients 
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relating transmission reconstruction and quantification and inferred anatomy-guided SPECT 
reconstruction and quantification. This increased correlation may be an artifact of the 
increased differences in the ROI means. Interestingly, the left frontal and central sulcus ROIs 
were also found to have marked rCBF deficits, suggesting that inferred anatomy may have 
difficulties near regions with substantially reduced radiopharmaceutical uptake. 

Figure 12 shows the means and standard errors for transmission reconstruction 
and inferred anatomy in each of the 12 ROIs, averaged across the entire population. The 
mean absolute difference for all ROIs across the whole population was 7.5%. Correlation for 
all ROIs and all patients was found to be high: r=0.92 as illustrated in Figure 13. 

TABLE 1 



Region 


ANOVA 


Paired sample 


Paired sample 


Significance 


correlation 


Significance 


Left frontal 


0.001 


0.965 


<0.001 


Right frontal 


0.088 


0.953 


<0.001 


Left central sulcus 


0.002 


0.957 


<0.001 


Right central sulcus 


0.234 


0.937 


<0.001 


Left parietal 


0.092 


0.772 


0.009 


Right parietal 


0.473 


0.862 


0.001 


Left temporal 


0.957 


0.939 


<0.001 


Right temporal 


0.397 


0.835 


0.003 


Left occipital 


0.271 


0.835 


0.003 


Right occipital 


0.073 


0.858 


0.001 


Left cerebellar 


0.545 


0.938 


<0.001 


Right cerebellar 


0.603 


0.925 


<0.001 



Error Propagation Analysis 

The results from the propagation analysis are presented in Table 2 below. 
These results are presented as errors relative to the "gold-standard" transmission 
reconstruction-guided reconstruction and spatial normaUzation. Three sources of error due to 
inferring anatomy were analyzed namely, errors in scatter distribution estimates, errors due to 
misguided attenuation compensation, and errors due to region-of-interest misplacement. The 
error components, averaged across ten patients, are shown for each of the twelve bilateral 
regions-of-interest. The table summarizes the percentage that each of these error sources 
contributed to the total quantitative differences between SPECT reconstructions guided by 
transmission scans or inferred anatomy. The fourth numerical column indicates the 
percentage of total differences that could not be accounted for by inferring anatomy. 
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On average, it was found that errors in scatter distribution estimates results in 
approximately 10.0% of the total quantitative error; attenuation correction: 36.6%, and ROI 
misplacement: 27.0%». The relative contributions of inferred anatomy-derived scatter and 
attenuation correction to the total error seem credible. Compensating for attenuation is of 
considerably greater consequence than removing scattered photons [Rosenthal et al 1995]. 
Approximately 26,5% of the total discrepancy between inferred anatomy and transmission 
imaging could not be accounted for in the error propagation analysis of scatter, attenuation, 
and ROI misplacement but is belived to be as a result of patient motion during data 
acquisition. 

TABLE 2 



Region 


Scatter 


Attenuation 


Region 


Unaccountable 


correction % 


correction % 


misplacement % 


error % 


Left frontal 


6.1 


44.5 


28.6 


20.8 


Right frontal 


12.0 


49.9 


29.4 


8.7 


Left central sulcus 


5.7 


40.4 


17.9 


36.0 


Right central sulcus 


6.3 


34.3 


16.0 


43.4 


Left parietal 


7.9 


25.7 


30.3 


36.1 


Right parietal 


10.0 


34.1 


25.9 


30.0 


Left temporal 


10.0 


38.3 


20.0 


3L8 


Right temporal 


13.8 


27.0 


30.7 


28.5 


Left occipital 


13.2 


23.2 


27.6 


36.0 


Right occipital 


8.3 


22.5 


29.0 


40.2 


Left cerebellar 


13.5 


48.7 


31.9 


5.8 


Right cerebellar 


12.5 


50.2 


36.4 


0.8 


Average 


10.0 


36.6 


27.0 


26.5 



Qualitative and Quantitative Comparisons 

Comparing inferred anatomy with transmission reconstructions indicated good 
reproduction of soft-tissue and hard-tissue in cortical areas for all ten patients. However, in 
many scans, a discrepancy was indicated near the sinus cavity, as shown in Figure 10. 
Despite this, inferred anatomy is more robust and accurate in providing estimates of 
underlying tissue distribution than fitting ellipses to photopeak emission data. The technique 
of fitting ellipses depends on adequate facial activity since the contours of the brain and head 
differ so considerably at the level of the cerebellum. Facial activity should not form the basis 
for estimating underlying tissue distributions since uptake varies with radiopharmaceutical 
and time between injection and scanning [Leveille et al 1992, Van Dyck et al 1996], making 
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it an unreliable dependency. Parenthetically, the qualitative similarities demonstrated 
between inferred anatomy and transmission reconstruction indicate confidence in accurately 
guiding scatter and attenuation correction. However, it is important to note that similar 
shape is a sufficient, but not necessary, prerequisite for accurate scatter and attenuation 
correction [Welch et al 1997, Natterer 1993]. 

Slight truncation effects are noticed on the transmission images for three of the 
ten scans. Truncation occurred for kyphotic patients with broad shoulders or with short 
necks. This limitation was generally restricted to transaxial slices below the level of the 
cerebellum or very near its base, where the gamma camera's circular field-of-view proved to 
be inconvenient. The truncation only involved nasal cartilage, and is therefore not expected 
to significantly impact results, as is demonstrated in the quantitative accuracy exhibited at the 
cerebellar level (see Figure 12). 

Inferred anatomy had difficulties in regions with marked rCBF deficit, such as 
the left frontal lobe. Although frontal lobes generally exhibit high variability in HMPAO 
uptake [Deutsch et al 1997], only the left frontal lobe had statistically significant error. This 
suggests that regional quantitative errors incurred by inferred anatomy are associated with 
rCBF deficits. However, previous work demonstrates that they may also be sensitive to 
spatial region. Achieving good quantitative accuracy in extended sources is often difficult 
and this is particularly true with peripheral ROIs, where reconstructed activity is most 
sensitive to misregistration of the attenuating medium [Huang et al 1979]. For example, in 
the brain a mismatch between emission and transmission data no greater than 5mm can 
produce a 10% error in a 10 mm thick peripheral cortical ROI [Andersson et al 1995b, Huang 
et al 1979]. In general, extended sources, such as the brain, are affected by misregistration 
more than compact sources, such as the heart [Andersson et al 1995b, McCord et al 1992, 
Ter-Pogossian 1992, Matsunari et al 1998]. Since the brain is elliptical, it is expected that 
regions along the periphery of the long axis of the head will be more sensitive to attenuation 
map misregistration than those along the short axis. 

Although a preferred embodiment of the present invention has been described, 
those of skill in the art will appreciate that variations and modifications may be made to the 
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present invention without departing from the spirit and scope thereof as defined by the 
appended claims. 
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APPENDIX A 

//initialize 3D array. Projections are 2D, but we have 
//many projections. 

//init 3D array. Store projections here after scatter 
//correction, 

//init 3D array. Store projections here after scatter and 
//attenuation correction. 

//init 3D array. The reconstructed 3D brain scan, 
//init 3D array. Transmission projections of anatomy, 
//init 3D array. The 3D distribution of anatomy goes 
//here . 

emis proj = do_patient_einission_scan; 

// first step is to acquire the patient emission scan. 
// The acquire data is in the form of projections of the 
// radiopharmaceutical distribution. 

// These projections require scatter and attenuation 
// correction. 

// In order to perform these corrections, an estimate of 
// anatomy must be obtained. Current, two methods of 
// obtaining an anatomy estimate are possible: 

if (infer==l) 

anatomy = inf er_anatoniy (emis_proj ) ; // This is our method. See the 

// function, below. This is what 
// we are patenting. 

elseif {infer==0) 

{ 

// Alternatively, we acquire transmission projections of 
// the patient. 

// This is very similar to x-ray CT. 
tx proj = do_patient_transmission_scan; 



einis_proj = 

emis_proj_sc = 

emispro j_sc_ac = 

emis_reco_sc ac = 
txproj = 
anatomy = 
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anatomy = recons truct ( tx__proj ) ; // These projections are 

// reconstructed into a 3D 

// distribution of anatomy. Note 

// that this reconstruction is 

// almost identical to the 

// emission reconstruction; 

// however, transmission 

// reconstructions do not 

// require scatter or attenuation 

// correction. 

} 

// After we have our emission projections and an anatomy 
// estimate, we sequentially apply scatter correction and 
// attenuation correction to the emission data. 

emispro j_sc = scatter^correction (emis_proj , anatomy) ; 

// Perform the scatter correction on the 
// emission_data . 

// The scatter correction requires knowing the 
// anatomy. 

emis_proj_sc_ac = attenuation_correction (emis_proj_sc , anatomy) ; 

// Perform the attenuation correction on the 
// emission data that has been scatter corrected, 
// The attenuation correction requires knowing the 
// anatomy. 

eiais_reco_sc_ac = reconstruct (emis_proj_sc_ac ) ; 

// After correcting the emission data for scatter 
// and attenuation, perform the emission 
// reconstruction. 

Stop 
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function aligned_head_atlas = inf er anatomy (emis_pro j ) ; 

{ 

emis^reco = // init 3D array. This will hold a preliminary emission 

// reconstruction that has NOT been corrected for scatter 
// or attenuation. 

head_atlas = // init 3D array. This will hold the full head 

// atlas (brain, skull, soft tissue) . The head 

// atlas is in its original orientation. 

head atlas_brain = // init 3D array. This will hold the brain 

// component of the head atlas. This brain 

// component will be in its original 

// orientation. 

alignment_par = //a set of parameters that represent the 3D 

// transformations (including, but not limited to: 
// rotation, shifting, scaling) necessary to align 
// head_atlas_brain to emis_reco. These 
// transformations will be applied to head_atlas . 
aligned_head_atlas = // init 3D array. The alignment parameters 

// are applied to head_atlas. This variable 
// is the function's output. 



emis_reco = reconstruct (eiais_proj ) // Perform a preliminary 

// reconstruction of the emission 
// data. 
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APPENDIX A 



head_atlas ~ load_3D_data (atlas_f ilename) ; // Load the head 

// atlas. 

head_atlas_brain = extract (head_atlas, brain) ; // Identify and 

// extract the brain 
// from the head 
// atlas. 

aligiunent_par = f ind_optimal_alignment (emis reco, head atlas brain) ; 

// Calculate the 3D transformation parameters 
// to align head_atlas__brain to emis_reco. 
// This procedure calculates the optimal set 
// of transformation parameters. In medical 
// science, this alignment is known as 
// "registration" . 

aligned_head_atlas = apply_alignment (head_atlas , alignment_par) ; 

// Apply the alignment to head_atlas . 
// aligned_head_atlas 

// represents an anatomy estimate that is 

// used for guiding the scatter and 

// attenuation corrections. (Note that the 

// variable aligned_head__atlas is passed back 

// as anatomy to the main program.) Also note 

// that, although the set of transformation 

// parameters that was calculated represents 

// the optimal set for aligning 

// head_atlas_brain to emis_reco, these same 

// parameters are now applied to head_atlas . 
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What is claimed is : 

I A method of applying scatter and attenuation correction to emission 

tomography images of a region of interest of a subject under observation comprising the steps 
of: 

aligning a three-dimensional computer model representing the density 
distribution within said region of interest with said emission tomography images; and 

applying scatter and attenuation correction to said emission tomography 
images using said aligned computer model as a guide. 

2. The method of claim 1 wherein said computer model is in the form of a two- 
component atlas. 

3. The method of claim 2 wherein during said aligning step, a functional 
component of said atlas is firstly aligned with said emission tomography images to generate a 
set of spatial transformation parameters and thereafter, the anatomical component of said atlas 
is aligned with said emission tomography images using said set of spatial transformation 
parameters. 

4. The method of claim 3 wherein said functional component simulates a SPECT 
or PET scan of said region of interest and wherein said anatomical component simulates a 
transmission scan of said region of interest. 

5. The method of claim 4 wherein said region of interest is the head and wherein 
said functional component is the brain component of a head atlas. 

6. The method of claim 4 wherein said region of interest is the heart, said 
functional component of said atlas simulating a cardiac image and said anatomical component 
of said atlas representing anatomical features of the thorax. 

7. The method of claim 6 wherein the anatomical features of the thorax include: 
soft-tissues such as the heart, liver, muscle, and fat; 
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very low-density soft-tissues such as the lungs; and 

high-density tissues such as bone and cartilage in the ribs and spine, 

8^ The method of claim 3 further comprising the step of selecting an atlas from a 

database of atlases prior to performing said aligning step. 

9 The method of claim 8 wherein said selecting step is performed manually. 

10. The method of claim 8 wherein said selecting step is performed automatically 
based on the degree of registration of each atlas in said database with said emission 
tomography images. 

1 1 . The method of claim 10 wherein the degree of registration is determined by: 
performing a preliminary reconstruction of each atlas; and 

registering the atlas to the preliminary reconstruction. 

12. The method of claim 10 further comprising the step of combining multiple 
atlases to yield a resuUant atlas that better registers with said emission tomography images. 

13 The method of claim 8 wherein said database includes disease specific atlases, 

physical trait specific atlases and/or tracer or lesion specific atlases. 

14. In an emission tomography imaging system where emission tomography 
images of a region of interest of a subject are taken for analysis and are corrected for scatter 
and attenuation, the improvement comprising: 

using a three-dimensional computer model approximating the density 
distribution within the region of interest as a guide to the application of scatter and 
attenuation correction. 

15. The emission tomography imaging system of claim 14 wherein said computer 
model is in the form of a two-component atlas. 
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16. The emission tomography imaging system of claim 15 wherein during said 
aUgning step, a functional component of said atlas is firstly aligned with said emission 
tomography images to generate a set of spatial transformation parameters and thereafter, the 
anatomical component of said atlas is aligned with said emission tomography images using 
said set of spatial transformation parameters. 

17. The emission tomography imaging system of claim 16 wherein said functional 
component simulates a SPECT or PET scan of said region of interest and wherein said 
anatomical component simulates a transmission scan of said region of interest. 

18. An emission tomography image processing system comprising: 
memory storing emission tomography images of a region of interest of a 

subject; 

said memory also storing at least one three-dimensional computer model of 
said region of interest, said computer model representing the density distribution within said 

region of interest; and 

a processor for registering said computer model with said emission 
tomography images and for applying scatter and attenuation correction to said emission 
tomography images using said registered computer model as a guide. 

19. An emission tomography image processing system as defined in claim 18 
wherein said at least one computer model is in the form of at least one two-component atlas. 

20. An emission tomography image processing system as defined in claim 19 
wherein said processor firstly registers a functional component of said atlas with said 
emission tomography images to generate a set of spatial transformation parameters and then 
registers the anatomical component of said atlas with said emission tomography images using 
said set of spatial transformation parameters. 

21. An emission tomography image processing system as defined in claim 20 
wherein said functional component simulates a SPECT or PET scan of said region of interest 
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and wherein said anatomical component simulates a transmission scan of said region of 
interest. 

22. An emission tomography image processing system as defined in claim 21 
wherein said memory stores a database of atlases. 

23. An emission tomography image processing system as defined in claim 22 
wherein said processor selects an atlas from said database automatically based on the degree 
of registration of each atlas in said database with said emission tomography images. 

24. An emission tomography image processing system as defined in claim 23 
wherein said processor performs a preliminary reconstruction of each atlas and registers the 
atlas to the preliminary reconstruction to determine the degree of registration of each atlas. 

25. An emission tomography image processing system as defined in claim 23 
wherein said processor combines multiple atlases to yield a resultant atlas that better registers 
with said emission tomography images. 

26. An emission tomography image processing system as defined in claim 22 
wherein said database includes disease specific atlases, physical trait specific atlases and/or 
tracer or lesion specific atlases. 

27. An emission tomography imaging system comprising: 

means for taking emission tomography images of a region of interest of a 
subject to form a three-dimensional image of said region of interest; 

memory to store said emission tomography images, said memory also storing 
at least one three-dimensional computer model of said region of interest, said computer model 
representing the density distribution within said region of interest; and 

a processor for aligning said computer model with said emission tomography 
images and for applying scatter and attenuation correction to said emission tomography 
images using said aligned computer model as a guide. 
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28. A computer readable medium including computer program code for applying 
scatter and attenuation correction to emission tomography images of a region of interest of a 
subject, said computer readable medium including: 

computer program code for aligning a three-dimensional computer model 
representing the density distribution within said region of interest with said emission 
tomography images; and 

computer program code for applying scatter and attenuation corrections to said 
emission tomography images using said aligned computer model as a guide. 

29. A computer readable medium as defined in claim 28 wherein said computer 
program code for ahgning includes: 

computer program code for aligning a functional component of said computer 
model simulating a SPECT or PET scan of said region of interest and for generating a set of 
spatial transformation parameters; and 

computer program code for aligning an anatomical component of said 
computer model simulating a transmission scan of said region of interest using said set of 
spatial transformation parameters. 
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SPECT scan axial slice 




Appurtenant anatomy for scatter 
and attenuation correction (sagittal) 
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SPECT scan without scatter 
and attenuation coiTcction 
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Transmission scan axial slice 




Appurtenant anatomy for scatter 
and attenuation correction (transverse) 



SUBSTITUTE SHEET (RULE 26) 



wo 00/10034 PCT/CA99/00751 



Figure 2 2/13 



50 



Patient 



60 



62 



Scatter and 
attenuation 
corrections 



Emission 
reconstruction 



64 



Patient scan 
analysis 



Emission 
reconstruction 



52 



i 



56 



Calculate alignment of brain component 
to emission reconstruction 



Alignment 
parameters 



1 



54 



Head 
atlas 



Apply alignment 
to entire atlas 



58 



SUBSTITUTE SHEET (RULE 26) 



wo 00/10034 




PCT/CA99/00751 



Figure 3 
Phantom data 




SPECT scan (transverse) Transmission scan (transverse) 




Appurtenant anatomy for 
scatter and attenuation correction 
(sagittal) 
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Appurtenant anatomy for 
scatter and attenuation correction 
(transverse) 





Inferred anatomy for 
scatter and attenuation correction 
(sagittal) 



Inferred anatomy for 
scatter and attenuation correction 
. (transverse) 
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Slice 80 Slice 94 

(A) Mid-cortical and (B) cerebellar transverse slices. 
Volume enclosed by brain contour (white), appurtenant 
anatomy (gray), and inferred anatomy (black). 
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Transverse slice 



Both the brain contour and uniform inferred anatomy differ 
from appurtenant anatomy. Here, this difference (in pixels) 
is plotted as a function of transverse slice. 
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Brain 
contour 




Sagittal slice through (A) appurtenant anatomy and 
(B) inferred anatomy. The (appurtenant) brain contour 
is shown as the thick black outline in both (A) and (B). 
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A 



SPECT scan of the head phantom, (transverse slice) 
Black squares in (A) are a subset of the two-hundred 
macrovoxels defined in the reconstructed volume. 
The macrovoxels are used in the error metrics 
defined below. 
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Global bias, shown in (B), was defined as 
[(mean macrovoxel activity/calibration)-l]*100%. 
Data is shown for appurtenant anatomy (AA), 
braincontour (BC), and inferred anatomy (lA). 
Uniform corrections are shown in gray, and 
non-uniform in black. 



AA BC lA 



Global uniformity, shown in (C), was defined as 
[(coefficient of variation of macrovoxels)* 100%]. 
Data is shown for appurtenant anatomy (AA), 
braincontour (BC), and inferred anatomy (lA). 
Uniform corrections are shown in gray, and 
non-uniform in black. 
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A 



(A) SPECT scan of phantom (transverse slice 
at cortical level). Hot and cold regions of interest 
were defined (hot at centre). The regions were 
used to calculate contrast recovery and 
relative quantification, as described 
in (B) and (C) below. 
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Contrast recovery (B) was calculated as 
[(Hot / Cold activity) / true activity ratio]* 100%. 
Data is shown for appurtenant anatomy (AA), 
brain contour (BC), and inferred anatomy (lA). 



AA BC lA 



Relative quantification (cortical: cerebellar 
activity ratio - regions of interest not shown) 
is presented in (C). Data is shown for 
appurtenant anatomy (AA), brain 
contour (BC),and inferred anatomy (lA). 
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The Inferred Anatomy algorithm requires a patient's preliminary 
reconstruction (A) and a head atlas (B). The preliminary reconstruction does 
not include scatter and attenuation correction. The head atlas consists of a 
functional component (brain), and an anatomical component (extra-cerebral 
tissue). The functional component is then spatially registered to the 
preliminary reconstruction (C), and this optimal transformation is recorded. 
The patient's anatomy is then inferred by applying the same transformation 
to the anatomical component of the head atlas (D). 
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C D 




Examples of a patient's SPECT scan (brain) superimposed onto the 
pertinent transmission reconstruction (A and C, extra-cerebral tissue) 
and Inferred Anatomy (B and D, extra-cerebral tissue). The 
reconstruction includes scatter and attenuation correction guided by 
either the transmission reconstruction or Inferred Anatomy. 
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100 




Position [cm] 



Profiles (B) through reconstructions at the cortical-axial level (A) 
comparing a SPECT reconstruction guided by the pertinent 
transmission scan (solid) with the reconstruction guided by Inferred 
Anatomy (dots). The profiles were taken between the two arrows in A. 
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120 




Regional ceebral blood flow with voxel-by-voxel Lassen's correction, relative 
to the cerebellum is shown for the twelve bilateral regions-of-interest. Means 
and standard errors are calculated across the ten patients. The quantification 
follows scatter and attenuation correction guided by transmission 
reconstructions (circles) or Inferred Anatomy (squares). 
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Comparison of reconstructions guided by transmission scans (abscissa), 
with reconstructions guided by Inferred Anatomy (ordinate). 120 data 
points are shown; these are from twelve regions-of-interest and ten 
patients. Each point represents regional cerebral blood flow with 
Lassen's correction, relative to the cerebellum. The equation of the line- 
of-best-fit is y=1.02x-3.83, and the correlation coefficient squared is 
r2=0.92 (p«0.0001). 
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The applicant is entitled, if he so wishes, to amend the claims of the International Application (see Rule 46): 

When? The time limit for filing such amendments is normally 2 months from the date of transmittal of the 
International Search Report; however, for more details, see the notes on the accompanying sheet. 

Where? Directly to the International Bureau of WlPO 
34, chemin des Colombettes 
121 1 Geneva 20. Switzerland 
Fascimile No.: (41-22) 740.14.35 

For more detailed instructions, see the notes on the accompanying sheet. 

2. I I The applicant is hereby notified that no Internationat Search Report will be established and that the declaration under 
' — ' Article 1 7(2)(a) to that effect is transmitted herewith. 

3. \^ With regard to the protest against payment of (an) additional fee(s) under Rule 40.2. the applicant is notified that: 

I I the protest together with the decision thereon has been transmitted to the International Bureau together with the 
' — ' applicant's request to forward the texts of both the protest and the decision thereon to the designated Offices. 

I I no decision has been made yet on the protest; the applicant will be notified as soon as a decision is made. 

4. Further action(s): The applicant is reminded of the following: 

Shortly after 18 months from the priority date, the international application will be published by the International Bureau. 
If the applicant wishes to avoid or postpone publication, a notice of withdrawal of the international application, or of the 
priority claim, must reach the International Bureau as provided in Rules SObisA and 90£j/s.3, respectively, before the 
completion of the technical preparations for internationat publication. 

Within 19 months from the priority date, a demand for international preliminary examination must be filed if the applicant 
wishes to postpone the entry into the national phase until 30 months from the priority date (in some Offices even later). 

Within 20 months from the priority date, the applicant must perform the prescribed acts for entry into the national phase 
before all designated Offices which have not been elected In the demand or in a later election within 19 months from the 
priority date or could not be elected because they are not bound by Chapter II. 



Name and mailing address of the International Searching Authority 
European Patent Office. P.B. 5818 Patentlaan 2 
y*l» NL-2280 HV Rijswijk 

S/JI Tel. (+31-70) 340-2040. Tx. 31 651 epo nl. 
Fax: (+31-70) 340-3016 



Authorized officer 

Emilio Fontana Balparda 



Form PCT/lSA/220 (July 1 998) 



INTERNATIONAL SEARCH REPORT 



A. CLASSIFICATION OF SUBJECTlWlfrER 

IPC 7 GOlTl/164 



^nter 



International Application No 

T/CA 99/00751 



According to International Patent Classification (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 7 GOIT 



Documentation searched other than minimum documentation to the ejrtent that such documents are included in the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category " Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



S.L. BACHARACH, ET AL: "Attenuation 

correction in cardi " 

JOURNAL OF NUCLEAR CARDIOLOGY, 

vol. 2, no. 3, June 1995 (1995-06), pages 

246-255, XP002125238 

St Louis, USA 

cited in the application 

page 247, column 2, line 5 -page 248, 

column 1, line 6; figure 2 



□ 



Further documents are listed in the continuation of box C. 



□ 



Patent family members are listed in annex. 



^ Special categories of cited documents : 

"A" document defining the general state of the art which is not 
considered to be of particular relevance 

"E" earlier document but published on or after the international 
filing date 

"L" document which may throw doubts on priority claim(s) or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

"O" document referring to an oral disclosure, use, exhibition or 
other means 

"P" document published prior to the intematlonal filing date but 
later than the priority date claimed 



T" later document published after the International filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory undertying the 
Invention 

"X° document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

"Y" document of particular relevance: the claimed invention 

cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 

"&" document member of the same patent family 



Date of the actual completion of the international search 



28 December 1999 



Date of mailing of the international search report 



13/01/2000 



Name and mailing address of the ISA 

European Patent Office, P.B. 5818 Patentlaan 2 
NL - 2280 HV Rijswijk 
Tel. (Hh31-70) 340-2040. Tx. 31 651 epo nl. 
Fax: (+31-70) 340-3016 



Authorized officer 



Anderson, A 



Form PCT/lSA/210 (second sheet) (July 1992) 



PATENT COOPERATION TRE 

PCT 



INTERNATIONAL SEARCH REPORT 

(PCT Article 18 and Rules 43 and 44) 



Applicant's or agent* s file reference 

10319-2 DAR 


POR FURTHER Notification of Transmittal of International Search Report 

(Form PCT/ISA/220) as well as, where applicable, item 5 below. 

ACTION 


International application No. 

PCT/CA 99/00751 


International filing date (day/month/year) 

13/08/1999 


(Earliest) Priority Date (day/month/year) 

14/08/1998 


Applicant 

PRATO, Frank, S. et al . 



This International Search Report has been prepared by this International Searching Authority and is transmitted to the applicant 
according to Article 18. A copy is being transmitted to the International Bureau. 



sheets. 



This International Search Report consists of a total of Z 

nn It is also accompanied by a copy of each prior art document cited in this report. 



1 . Basis of the report 

a. With regard to the language, the international search was carried out on the basis of the international application in the 
language In which it was filed, unless otherwise indicated under this item. 

I I the international search was carried out on the basis of a translation of the international application furnished to this 
Authority (Rule 23-1 (b)). 

b. With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the international search 
was carried out on the basis of the sequence listing : 

I I contained in the international application in written form. 

filed together with the international application in computer readable form, 
furnished subsequently to this Authority in written form, 
furnished subsequently to this Authority in computer readble form. 



2. 
3. 



□ 
□ 
□ 
□ 

□ 

□ 
□ 



the statement that the subsequently furnished written sequence listing does not go beyond the disclosure in the 
international application as filed has been furnished. 

the statement that the information recorded in computer readable form is identical to the written sequence listing has been 
furnished 

Certain claims were found unsearchable (See Box t). 
Unity of invention is lacking (see Box It). 



4. With regard to the title, 

PC] the text is approved as submitted by the applicant. 

I I the text has been established by this Authority to read as follows: 



With regard to the abstract, 

PC] the text is approved as submitted by the applicant. 

I I the text has been established, according to Rule 38.2(b). by this Authority as it appears in Box III. The applicant may, 
' — ' within one month from the date of mailing of this international search report, submit comments to this Authority. 

The figure of the drawings to be published with the abstract is Figure No. 2 



Pn as suggested by the applicant. Q None of the figures. 

I I because the applicant failed to suggest a figure, 

I I because this figure better characterizes the invention. 



Form PCT/lSA/210 (first sheet) (July 1998) 



# 



NOTES TO FORM PCT/ISA/220 (c nt 



The letter must indicate the differences between the claims as filed and the claims as amended. II must, in 
particular, incftcate, in connection with each claim appearing in the international application (it being understood 
that identical indications concerning several claims may be grouped),whether 



(0 


the claim 


ia unchanged; 


00 


the claim 


is cancelled; 


Oii) 


the claim 


ts new; 


Ov) 


the claim 


replaces one or more daims as filed; 


(V) 


the claim 


is the result of the division of a claim as filed 



The following examplas Illustrate the manner In which amendments must be explained In the 
accompanying letter: 

1 . [Where originally there were 48 claims and after amendment of some claims there are 51 ]: 
'Qaims 1 to 29, 31 , 32. 34, 35, 37 to 48 replaced by amended claims bearing the same numbers; 
daims 30, 33 and 36 unchanged; new daims 49 to 51 added." 

2. [Where originally there were 15 claims and after amendmerrt of all daims there are 11]: 
'Claims 1 to 1 5 reptaoed by amended daims 1 to 1 1 

3. [Where originally there were 1 4 claims and the amendments consist in cancelling some daims and in ad<£ng 
r>ew claims]: 

"Claims 1 to 6 and 14 unchanged; claims 7 to 13 cancelled; new claims 15, 16 and 17 added.' or 
'Claims 7 to 13 cancelled; new claims 15, 16 and 17 added; all other claims unchanged.* 

4. [Where various kinds of amendments are made]: 

'Claims 1 -1 0 unchanged; claims 1 1 to 1 3, 1 8 and 1 9 cancelled; claims 1 4, 1 5 and 1 6 replaced by amended 
daim 14; claim 17 subdivided into amertded daims 15, 16 and 17; new daims 20 and 21 added.* 



"Statement under article 19(1)** (Rule 46.4) 

The amendments may be accompanied by a statement explaining the amendments and indcating any impact 
that such amer>dments might have on the description and the drawings (which cannot be amended under 
Artide19(1)). 

The statement will be published virith the irrtemational application and the amended daims. 
It must be In the language In which the International apppHcaticn Is to be published. 

It must be tsrief, not exceeding 500 words if in English or if translated into English. 

It should not be confused with and does not replace the letter indicating the cfifferenoes between the claims 
as filed and as amended. It must be filed on a separate sheet and must be identified as such t>y a heading, 
preteratsly by using the words "Statement under Article 19(1).* 

It may not contain any disparaging commerrts on the international search report or the relevance of citations 
contained in that report. Reference to citations, relevant to a given daim, contairted in the international search 
report may be made only in connection with an amendment of that claim. 



Consequence if a demand for International preliminary examination has already been filed 

If, at the time of filing any amendments under Article 1 9, a demand for international preliminary examination 
has already been submitted, the applicant must preferat^ty, at the same time of filing the amendments with the 
International Bureau, also file a copy of such amendments with the Intemational Preliminary Examining 
Authority (see Rule 62.2(a), first sentence). 



Consequence with regard to translation of the tntematlonal application for entry Into the national phase 

The applicant's attention is drawn to the fact that, where upon entry into the national phase, a trarwlation of the 
daims as amended under Article 19 may have to be furnished to the designated/elected Offices, instead of. or 
in addition to, the translation of the daims as filed. 

For further details on the requirements of each designated/elected Office, see Volume II of the PCT Applicant's 
Guide. 



Notes to Form PCT/ISA/220 (second sheet) (January 1994) 
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I^I^S TO FORM PCT/ISA/220 



These Notes are intended to grve the basic instructions concerning the filing of amendments under article 19. The 
Notes are based on the requirements of the Patent Cooperation Treaty, the Regulations and the Administrative Instmctiorw 
under that Treaty. In case of discrepancy between these Notes and those requirements, the latter are applicable. For more 
detailed information, see also the PCT Applicant's Guide, a publication of WtPO. 

!n these Notes, 'Article', "Rule", and 'Section' refer to the provisions of the PCT, the PCT Regulations and the PCT 
Administrative Instructions respectively. 



INSTRUCTIONS CONCERNING AMENDMENTS UNDER ARTICLE 19 



The applicant has, after having received the intemational search report, one opportunity to amend the claims of the 
intemationaJ application. It should however be emphasized that, since all parts of the intemational appticaiion (claims, 
description and drawings) may be amended during the intemational preliminary examination procedure, there is usually 
no need to file amendments of the claims under Article 19 except where, e.g. the applicant wants the latter to be published 
for the purposes of provisional protection or has another reason for amending the claims before international pbulication. 
Furthermore, it should t>e emphasized that provisional protection is available in some States only. 



What parts of ttie Intemationai applicadon may be amended? 

Urrder Article 1 9, only the claims may be amended. 

During the intemational phase, the daims may also fc^e amended (or further amended) under Article 34 before 
the International Preliminary Examining Authority. The description and drawings may only be amended under 
Article 34 before the International Examining Authority. 

Upon entry into the national phase, alt parts of the intemationaJ application may be amended under Article 28 
or, where applicable, Article 41 . 



When? Within 2 months from the date of trarrsmittal of the intemational search report or 1 6 months from the priority 

date, whichever time limit expires later. It should be noted, however, that the amendments will be considered 
as having been received on time if they are received by the Intemational Bureau after the expiration of the 
applicable time limit but before the comptetion of the technical preparations for intemational publication 
(Rule 46.1). 



Where not to file the amendments? 

The amendmerrts may only be filed with the International Bureau and not with the receiving Office or the 
International Searching Authority (Rule 46.2). 

Where a demand for intemational preliminary examination haa been /is filed, see below. 



How? Either by canceHing one or more entire claims, by adding one or more new claims or by amending the text of 

one or more of the daims as filed. 

A replacement sheet must be submitted for each sheet of the dainns which, on account of an amendment or 
amendments, differs from the sheet originally filed. 

All the daims appearing on a replacement sheet must be numbered in Arabic numerals. Where a claim is 
cancelled, no renumbering of the other claims is required. In alt cases where dainrs are renumbered, they must 
k» renumbered consecutively (Administrative Instructions, Section 205(b)}. 

The amendments must be made in the language In which the Intemational application is to be published. 



What documents must/may accompany ttie amendments? 
Letter (Section 205(b)): 

The amendments must be submitted with a letter. 

The letter will not be published with the intemational application and the amended daims. It should not be 
confused with the "Statement under Article 1 9(1)* (see below, under "Statement under Article 19(1)*). 

The letter must be In English or French, at the choice of the applicant. However, If the language of the 
intemational application Is English, the letter must be In English; If the language of the Intemational application 
Is French, the tetter must be In French. 



Notes to Form PCT/1SA;220 (first sheet) (January 1 994) 
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PATENT CCX)PERAT10N TREATY 



From the 

INTERNATIONAL PRELIMINARY EXAMINING AUTHORITY 



To: 

RUSTON, D. 
Sim & McBumey 
330 University Avenue 
6th floor 

Toronto, Ontario MSG 1 R7 
CANADA 



ECi 



PCT 



NOTIFICATION OF TRANSMITTAL OF 
THE INTERNATIONAL PRELIMINARY 
EXAMINATION REPORT 

(PCT Rule 71.1) 



Date of mailing 
(day/month/year) 



13.09.2000 



Applicanfs or agent's file reference 
10319-2 DAR 


IMPORTANT NOTTFICATTON 


International application No. 
PCT/CA99/00751 


International filing date (day/month/year) 
13/08/1999 


Priority date (day/month/year) 
14/08/1998 


Applicant 

PRATO, Frank, S. et al. 



1 . The applicant is hereby notified that this International Prelinninary Exannining Authority transmits herewith the 
international preliminary examination report and its annexes, if any, established on the international application. 



2. A copy of the report and its annexes, if any, is being transmitted to the International Bureau for communication 
to all the elected Offices. 



3. Where required by any of the elected Offices, the Intemational Bureau will prepare an English translation of the 
report (but not of any annexes) and will transmit such translation to those Offices. 



4. REMINDER 

The applicant must enter the national phase before each elected Office by performing certain acts (filing 
translations and paying national fees) within 30 months from the priority date (or later in some Offices) (Article 
39(1)) (see also the reminder sent by the Intemational Bureau with Fomn PCT/IB/301). 

Where a translation of the intemational application must be furnished to an elected Office, that translation must 
contain a translation of any annexes to the intemational preliminary examination report. It is the applicant's 
responsibility to prepare and f umish such translation directly to each elected Office concerned. 

For further details on the applicable time limits and requirements of the elected Offices, see Volume II of the 
PCT Applicant's Guide. 



Name and mailing address of the IPEA/ 



European Patent Office 
D-80298 Munich 



Tel. +49 89 2399 - 0 Tx: 523656 epmu d 
Fax: +49 89 2399 - 4465 



Authorized officer 
DEL FRATE, A 

Tel.+49 89 2399-2413 




Form PCT/IPEA/416 (July 1992) 



ATENT COOPERATION T\ 




PCT 



INTERNATIONAL PRELIMINARY EXAMINATION REPORT 

(PCT Article 36 and Rule 70) 



Applicants or agents file reference 


See Notification of Transmittal of International 


10319-2 DAR 


FOR FURTHER ACTION preliminary Examination Report (Form PCT/IPEA/416) 


Internationa) application No. 


International filing date (day/month/year) 


Priority date (day/month/year) 


PCT/CA99/00751 


13/08/1999 


14/08/1998 


International Patent Classification (IPC) or national classification and IPC 




G01T1/164 






Applicant 






PRATO, Frank, S. et al. 







1 . This international preliminary examination report has been prepared by this International Preliminary Examining Authority 
and is transmitted to the applicant according to Article 36. 



2. This REPORT consists of a total of 4 sheets, including this cover sheet. 

S This report is also accompanied by ANNEXES, i.e. sheets of the description, claims and/or drawings which hav 
been amended and are the basis for this report and/or sheets containing rectifications made before this Authority 
(see Rule 70.1 6 and Section 607 of the Administrative Instructions under the PCT). 

These annexes consist of a total of 3 sheets. 



3. This report contains indications relating to the following items: 



II 


□ 


III 


□ 


IV 


□ 


V 




VI 


□ 


VII 




VIII 





Basis of the report 



citations and explanations suporting such statement 



Date of submission of the demand 
13/03/2000 



Date of completion of this report 



13.09.2000 



Name and mailing address of the international 
preliminary examining authority: 
European Patent Office 

D-80298 Munich 
Tel. +49 89 2399 - 0 Tx: 523656 epmu d 

Fax: +49 89 2399 - 4465 



Authorized officer 
Rabenstein, W 

Telephone No. +49 89 2399 2450 




Form PCT/IPEA/409 (cover sheet) (January 1994) 



INTERNATIONAL PRELIMINARY 
EXAMINATION REPORT 



International application No. PCT/C A99/0075 1 



1. Basis tth report 

1 . This report has been drawn on the basis of {substitute sheets which have been furnished to the receiving Office in 
response to an invitation under Article 14 are referred to in this report as "originally filed" and are not annexed to 
the report since they do not contain amendments,): 

Description, pages: 

1 ,2,4-30 as originally filed 

3 as received on 22/08/2000 with letter of 14/08/2000 



Claims, No.: 

1-6,7 (part),21 (part), as originally filed 
22-29 



7 (part) ,8-20 



as received on 



22/08/2000 with letter of 



14/08/2000 



Drawings, sheets: 

1 /1 3-1 3/1 3 as originally filed 



2. The amendments have resulted in the cancellation of: 

□ the description, pages: 

□ the claims, Nos.: 

□ the drawings, sheets: 

3. □ This report has been established as if (some of) the amendments had not been made, since they have been 

considered to go beyond the disclosure as filed (Rule 70.2(c)): 



4. Additional observations, if necessary: 



Form PCT/lPEA/409 (Boxes l-VIII, Sheet 1) (January 1994) 



INTERNATIONAL PRELIMINARY 
EXAMINATION REPORT 



International application No. PCT/CA99/00751 



V. Reasoned statement under Artici 35(2) with regard to novelty, inv ntiv step or industrial 
applicability; citations and explanations supporting such statement 

1. Statement 



Novelty (N) 


Yes: 


Claims 


1-29 




No: 


Claims 




Inventive step (IS) 


Yes: 


Claims 


1-29 




No: 


Claims 




Industrial applicability (lA) 


Yes: 


Claims 


1-29 




No: 


Claims 





2. Citations and explanations 
see separate sheet 

VII. Certain defects in the international application 

The following defects in the form or contents of the intemational application have been noted: 
see separate sheet 

VIII. Certain observations on the international application 

The following observations on the clarity of the claims, description, and drawings or on the question whether the 
claims are fully supported by the description, are made: 

see separate sheet 



Form PCT/lPEA/409 (Boxes I- VIII, Sheet 2) (January 1994) 



INTERNATIONAL PRELIMINARY International application No. PCT/CA99/00751 
EXAMINATION REPORT - SEPARATE SHEET 



1 Re It mV 

1.1 The following document (D) is referred to in this communication: 

D1: S.L. BACHARACH, ET AL: 'Attenuation correction in card!' JOURNAL OF 
NUCLEAR CARDIOLOGY, vol. 2, no. 3, June 1995 (1995-06), pages 246- 
255, XP0021 25238 St Louis, USA. 

1 .2 The application relates to emission tomography and a method of applying attenu- 
ation correction to emission tomography images. In the prior art transmission 
images are used to correct the emission images (D1). The invention is based on 
the finding that using a three dimensional computer model for the correction leads 
to similar results. Thereby, the obvious disadvantages of the prior art method, 
which required an additional image to be taken, are avoided. As the available prior 
art relates to the use of transmission images for attenuation correction, the subject 
matter of independent claims 1, 14. 18 and 27 is novel and involves an inventive 
step, so that the requirements of Art. 33(2) and (3) PCT are met. 

1 .3 The other claims depend on claims 1 , 1 4, 1 8 and 27 and therefore fulfill the 
requirements as well. 

2 Re Item Vlll 

2.1 Although claims 1 and 14 (18 and 27) have been drafted as separate independent 
claims, they appear to relate effectively to the same subject-matter and to differ 
from each other only with regard to the definition of the subject-matter for which 
protection is sought and/or in respect of the terminology used for the features of 
that subject-matter. The aforementioned claims therefore lack conciseness. 
Hence, claims 1. 14, 18 and 27 do not meet the requirements of Article 6 PCT. 

2.2 Claim 16 refers back to claims 15, 14 but mentions an aligning step which is not 
defined in these claims. Therefore, claim 16 is unclear (Art. 6 PCT). 

2.3 By replacing two sheets of the claims pages, one feature of claim 7 is repeated 
twice and the first part of claim 21 is lacking. 



Form PCT/Separate Sheet/409 (Sheet 1) (EPO-April 1997) 
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The atlas may be selected from a database of atlases based on degree of 
registration with the emission tomogr^hy images. Alternatively, multiple atlases maybe 
combined to yield a resultant atlas which better registers with the emission tomography 
images. 

According to another aspect of the present invention there is provided an 
emission tomogr^hy imaging method where emission tomography images of a region of 
interest of a subject are taken for analysis and are corrected for scatter and attenuation, the 
method further comprising the step of: 

using a three-dimensional computer model approximating the density 
distribution within the region of interest as a guide to the appUcation of scatter and 
attenuation correction. 

According to stiU yet another aspect of the present invention there is provided 
an emission tomogr^hy image processing system comprising: 

memory storing emission tomography images of a region of interest of a 

subject; 

said memory also storing at least one three-dimensional computer model of 
said region of interest, said computer model representing the density distribution within said 
region of interest; and 

a processor for registering said computer model with said emission 
tomography images and for applying scatter and attenuation correction to said emission 
tomography images using said registered computer model as a guide. 

According to still yet another aspect of the present invention there is provided 
an emission tomography imaging system comprising: 

means for taking entussion tomography images of a region of interest of a 
subject to form a three-dimensional image of said region of interest; 

memory to store said emission tomography images, said memory also storing 
at least one three-dimensional computer model of said region of interest, said computer model 
representing the density distribution within said region of interest; and 
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soft-tissues such as the heart, liver, muscle, and fat; 
very low-density soft-tissues such as the lungs; and 
high-density tissues such as bone and cartilage in the ribs and spine. 

8. The method of claim 3 further comprising the step of selecting an atlas from a 
database of atlases prior to performing said aligning step. 

9. The method of claim 8 wherein said selecting step is performed manually. 

10. The method of claim 8 wherein said selecting step is performed automatically 
based on the degree of registration of each atlas in said database with said emission 
tomography images. 

1 1 . The method of claim 10 wherein the degree of registration is determined by: 
performing a preliminary reconstmction of each atlas; and 

registering the atlas to the preliminary reconstruction. 

12. The method of claim 10 further comprising the step of combining multiple • 
atlases to yield a resultant atlas that better registers with said emission tomography images. 

13. The method of claim 8 wherein said database includes disease specific atlases, 
physical trait specific atlases and/or tracer or lesion specific atlases. 

14. An emission tomography imaging method where emission tomography images 
of a region of interest of a subject are taken for analysis and are corrected for scatter and 
attCTuation, the method further comprising the step of: 

using a three-dimensional computer model approximating the density 
distribution within the region of interest as a guide to the application of scatter and 
attenuation correction. 
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15, The emission tomography imaging method of claim 14 wherein said computer 

model is in the form of a two-component atlas. 



16. The emission tomography imaging method of claim 15 wherein during said 
aligning step, a functional component of said atlas is firstly aligned with said emission 
tomography images to generate a set of spatial transformation parameters and thereafter, ttie 
anatomical component of said atlas is aligned with said emission tomography images using 
said set of spatial transformation parameters. 

17. The emission tomography imaging method of claim 16 wherein said 
functional component simulates a SPECT or PET scan of said region of interest and wherein 
said anatomical component simulates a transmission scan of said region of interest. 

18. An emission tomography image processing system comprisiag: 
memory storing emission tomography images of a region of interest of a 

subject; 

said memory also storing at least one three-dimrasional computer model of 
said region of interest, said computer model representing the density distribution within said 
region of interest; and 

a processor for registering said computer model with said emission 
tomography images and for applying scatter and attenuation correction to said emission 
tomography images using said registered computer model as a guide. 

1 9. An emission tomography image processing system as defined in claim 1 8 
wherein said at least one computer model is in the form of at least one two-component atlas. 



20. An emission tomography image processing system as defined in claim 1 9 

wherein said processor firstly registers a functional component of said atlas with said 
emission tomography images to generate a set of spatial transformation parameters and then 
registers the anatomical component of said atlas with said emission tomography images using 
said set of spatial transformation parameters. 
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